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Abstract 

We calculate the magnetic polarization tensor of the photon and of the gluons in 
a two flavor color superconductor with a LOFF pairing that consists of a single 
plane wave. We show that at zero temperature and within the range of the values 
of the Fermi sea displacement that favors the LOFF state, all the eigenvalues of 
the magnetic polarization tensor are non-negative. Therefore the chromomagnetic 
instabilities pertaining to a gapless color superconductor disappear. 
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1 Introduction 



The LOFF ( Larkin-Ovchinnikov-Fudde-Ferrell ) state [1], [2] was introduced as a can- 
didate superconducting state of an alloy containing paramagnetic impurities with a fer- 
romagnetic spin alignment. The presence of a nonzero expectation value for the spins of 
the impurities and its corresponding coupling to the electrons gives rise to a displacement 
of the Fermi surfaces of the pairing electrons since they have opposite spins. The char- 
acteristic feature of this state is that the corresponding order parameter is not constant 
but depends on the coordinates. Though the experimental evidence for the LOFF state 
is still inconclusive, the subject has experienced a revived interest due to its relevance to 
the physics of color superconductors and compact stars [3]. 

The interaction between two quarks that renders quark matter at ultra high baryon 
density unstable to the formation of Cooper pairs and leads to the emergence of a color 
superconducting state, results from the perturbative one-gluon exchange in the color anti- 
symmetric diquark channel. All quark flavors, u, d and s can be treated as massless since 
the Fermi energies of the quarks are much higher than Aqcd- In contrast at moderately 
high baryon densities, say the density inside a compact star, where the Fermi energies 
of the quarks are only few hundred MeV, instanton effects are expected to dominate the 
interaction responsible for the pairing of quarks [4] , again in the color antisymmetric chan- 
nel. Since the main pairing channel is between quarks of different flavors, factors such as 
the large mass of the s quark and the different magnitudes of electric charges of the u 
and d quarks together with the color /electric neutrality condition induce a considerable 
displacement between the Fermi momenta of the pairing quarks. This makes the search 
for a realistic color superconductivity phase challenging. 

The LOFF state is one of several states that were proposed as a possible ground 
state for quark matter at moderately high baryon density [5], [6]. Others include the 
homogeneous gapless states (g2SC or gCFL ) [7], [8], [9] and the heterogeneous mixture 
of BCS states without Fermi sea displacement and normal phases [12]. The gapless state 
is the analog of the Sarma state [10] for quark matter and a potential ground state since 
the Sarma instability can be removed by imposing the charge neutrality constraint [7], 
[11]. But a subsequent calculation revealed that the squares of the Meissner masses of 
the gluons were all negative signaling a chromomagnetic instability in the presence of the 
gauge field [13], [14], [15], [16]. This instability represents a major problem that needs to 
be resolved before we would be able to identify the ground state of the quark matter. 

The mixed state on the other hand is free from chromomagnetic instabilities but it 
has the drawback that it is difficult to make a quantitative comparison of its free energy 
with one of the homogeneous CSC state. 

In a previous paper [15], we explored the relationship between the chromomagnetic 
instability and the LOFF state. For a g2SC, we showed explicitly that the chromomagnetic 
instability associated with the 8th gluon implies the existence of a LOFF state with a 
lower free energy and a positive Meissner mass square for the 8th gluon excitation. But 
the possibility remains that the chromomagnetic instability will emerge when we calculate 
the Meissner masses of the gluons that are associated with the 4-7th generators of the 
color group. In this paper we shall continue this project by calculating the Meissner 
masses of all the gluons of a LOFF state. Our result shows that the chromomagnetic 
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instability is completely removed at zero temperature and in the LOFF window where 
the LOFF state is energetically favored for a given Fermi momentum displacement. 

This paper is organized as follows: The free energy function of a two flavor LOFF CSC 
will be derived in section II. A discussion of its properties will also be included. In section 
III, the Meissner masses of all the gluonic excitations will be calculated. The analytic 
proof of the absence of the chromomagnetic instabilities for a small gap parameter A, at 
T = will be presented in section IV and the numerical results for an arbitrary A will be 
discussed in section V, where we extend the chromagnetic instability free region to the 
entire LOFF window at T = 0. Finally, in section VI will summarize our results. Some 
technical details will be presented in the appendices. Unless explicitly stated, we shall 
work exclusively with the single plane wave ansatz for which analytical investigation can 
be carried out quite far. We shall also designate the symbol 2SC-LOFF for a two flavor 
LOFF color superconducting state. 



2 The Free Energy Function of a 2SC-LOFF 

We consider the pairing of two quark flavors, u and d. The NJL effective Lagrangian in 
the massless limit is given by [17], 

+ G D (4>cl5e c T2ip)(4>l5e c T 2 ipc) (1) 

where xjj represents the quark fields and ipc — Cip is its charge conjugate with C = 17274. 
All gamma matrices are hermitian, ^ c ) mn = e cmn is a 3 x 3 matrix acting on the red, 
green and blue color indices and the Pauli matrices f act on the u and d flavor (isospin) 
indices. The chemical potential is also a matrix in color-flavor space, i.e. 

fi = fi-8T3 + S'Xs (2) 

with \[ being the Zth Gell-Mann matrix. 

The most general expression of the order parameter of 2SC-LOFF is a superposition 
of plane waves of diquark condensate, i.e. 

< ^(r)75A 2 r 2 ^ c (r) >= £ $<?e 2 ^, (3) 

In this article, we shall be mainly working with the type of the LOFF state that contains 
only one plane wave, more specifically the form of the order parameter will be given by 

< ^75A 2 T 2 ^c(r) >= Qe 2 ** (4) 

The condensate eq. (4) is invariant under a translation accompanied by a phase shift. 
Since the free energy depends on the magnitude of the condensate only, it is translation- 
ally invariant. But this is not clearly the case with the general form of the 2SC-LOFF 
condensate eq. (3). 
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Let's introduce a transformed quark field, 

X(r) = e-^>(f). (5) 
The NJL Lagrangian in terms of the new quark field reads 

£ = -xi^-r^r ■ qx + xiA^x + G s [(xx) 2 + (x^x) 2 ] 

+ G D (xcl5e c T 2 x)(xi^ c T2Xc), (6) 

while the condensate (4) takes the form 

< x(^)75A 2 r 2 xc(r) >= (7) 

By ignoring the chiral condensate and expanding the NJL Lagrangian to linear order 
in the fluctuation 

^75A 2 r 2 V ; c - (8) 
we derive the following mean field expression for the NJL Lagrangian 

A 2 dx 1 

CmF = ~4G^ ~ X1 ^dx~ ~ ' ^ X + X74/iX + 2 A ( _Xc75A2T2X + X7sA 2 t 2 xc), (9) 

where the gap parameter is given by A = 2Gd& > 0. The free energy T is given by the 
Euclidean path integral 



jr = _k BT i n [ J [d^pdij] exp ( J #xCmf)\ (10) 



with < X4 < (/cbT)^ 1 and will be minimized at equilibrium. We find it convenient 
to work with the difference of the free energy density between the superconducting and 
normal phases in a grand canonical ensemble at the same 5 and 5', 

r = ^-^|A=o), (ii) 

where Q is the volume of the system. The equilibrium conditions read 

dT dT 

= 0, -7T- = 0, (12) 



and 



It follows then that 



at the minimum. 



d 2 T 

> 0, 



<9A 2 ' dq 

a 2 r d 2 r 



(<9A 2 ) 



(8A 2 ) 2 8A 2 dq 

a 2 v a 2 r 

dA 2 dq dq 2 

d 2 v 



> 0. (13) 



dq 2 



> (14) 
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In terms of the Nambu-Gorkov representation of the quark field, 



''' ( T * Y . ) * = (\-\r») (15) 



and its Fourier transform, 



W) = v^E / ( 16 ) 

we find that 



J ACmf = n[ - _ + 5 ? / ^*p5- ( p)*p] (it) 

where iS _1 (P) is the momentum representation of the inverse quark propagator 



S~\P) 



f - rf ■ q + /i7 4 - 5^ a t 3 + <5'7 4 A 8 A75A 
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y -A75A2 f + tf-q- /i74 - <57 4 r 3 - 5'7 4 A 8 J 

= f - <57 4 r 3 + (-^7 • g + /i7 4 + (5 / 7 4 A 8 )p 3 + iA7 5 A 2 p 2 - (18) 



P = — r^P^, P = (p, —v) with v = (2n + l^/c^T being the Matsubara frequency, and 
we have introduced the Pauli matrices with respect to NG blocks, i.e. 

Pl= (°i 0) p2= {u ~o) p3= {o -/)• (19) 

The product of Dirac, color, flavor and NG matrices is a direct product and the total 
dimensionality of S _1 (P) is 4(Dirac) x 3(color) x 2(flavor) x 2(NG) = 48. 
Carrying out the path integration in Eq. (10), we find 

r-^-^y-/^i n det ^ p ) (20) 

AG D 2 (27r)3 ln det5- 1 (P)' [ ) 

where S(P) = <S(P)|a=o is the free quark propagator. 

The typical energy scale where the pairing occurs is characterized by the zero tem- 
perature gap parameter A in the absence of a displacement, i.e. 5 = 5' = 0. At weak 
coupling, A << ft and the pairing interaction extends over a width of u on both sides 
of the Fermi level with A « lu « fi. Furthermore the displacement 5' that is induced 
by the condensate, is of higher order and can be neglected. As a result we have the only 
displacement parameter 5 ~ Ao- Under this approximation, the quark propagator takes 
the approximate form 

c( P \ „ (-tf-P + Psld [\P- Psq\ - p. + Psjiv - St 3 ) + AA 2 7s7 4 ] , 

1 2 ( p - jl )2 + ( iu _ 5r3 _p.$ ) 2 + A 2 y ) 

while the determinant in eq. (20) is given by 

detS-^P) = (2p) 24 nt(P - P? ~ (™ + s5-p- q) 2 + (AA) 2 ] 2 (22) 

s,A 
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where s — ±1 and A = ±1,0 refer to the eigenvalues of r 3 and A 2 respectively. A more 
detailed derivation of eqs. (21) and (22) can be found in appendix A. It follows from the 
analytic continuation of (21) to the real energy iv — > p , that the low- lying excitation 
spectrum reads 

E$ = y/(p-ji)* + A*±6-p-q (23) 

for red and green colors and E$ — \p — ft\ ± 5 — p ■ q for blue color. The excitations (23) 
are fully gapped if A > 5 + q. E-~ becomes gapless when \5 — q\ < A < 5 + q and both 
Ep become gapless when q > 5 and q — 5 < A < q + § [2] . 
Substituting (22) into (20) we obtain that 

^ A 2 p 2 k B T [i too e-z 2 (v,0,A) 

^ dcos9 d£ln (24) 



AG D vr 2 ^J-i J-oo s ?-z*(v,6,0) 



with £ = p — p, and z(v, 9, A) = y (iv + 5 — q cos 9) 2 — A 2 . The contributions from the 
eigenvalues of T3 with opposite sign are equal if we reverse the signs of the Matsubara 
frequency and simultaneously perform the transformation 9 — > n — 9. This will also be 
the case with the calculation of the Meissner masses in the next section. 
We shall employ the phase conventions that are determined by 

lmz(u, 9, A) > for v > , (25) 

with ((—u,9, A) = (*(u,9, A). Let's now proceed and perform the £ integral by closing 
the contour around one of the branch cuts of the logarithm. As a result we obtain the 
expression 

r.2l 



r = -^-^^Im f dcos9 £ [z(u,9,A)-z(u,9,0)]- (26) 
At zero temperature, T = 0, we have 

« B T £ (...) = £•£(...) (27) 



0<U<OJQ 



and 

7T 2 V 



A 2 2/i 2 / 2cj 1 



AG D n 2 V A 2 



+ rfcos4g(|^|-A)(A 2 ln ' 91 Y |W5 fl 2 -A 2 ) + 5fl (28) 

TT Z J-l 1 A v J 

where 5g = (1 — pcos9)5 with p — q/8, being the normalized diquark momentum. The 
limit of the angular integral depends on the gapless region of the excitation spectrum (23) 
and will be discussed in detail in the appendix B. In terms of the gap parameter A at 
5 = 0, given by the standard BCS gap equation 

1 2p2 m^, (29) 



4G D 7T 2 A ' 
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we have 



A 



2 ) + 



(P + 1) 3 5 2 



4p 



\±^ + \(2x\-?>x l + l) 
1 — xi 3 



+ 



(P~l) ; 
4p 



(1 - x\) In + ^ - 3x 2 + 1)1 }, (30) 

1 — X 2 O J 



where Xi and x 2 are the same dimensionless parameters introduced in Ref. [2], i.e. 

A 



and 



Xi 



%2 



9(1- 



1 



A 2 



(p+l)6'\ (p+l) 2 <5 2 
A 



p-l\6'\ (p-l) 2 S r 



(31) 



(32) 



Differentiating eq. (30) with respect to A 2 and p, we derive the equilibrium equations 



or 

M 2 



^ 2 r ln A 

A 



7T 



4p 



X\ 



1 +Xi 



2xi 



P-1/, 1-X 2 



4p 



In 



l + x 2 



2x< 



(33) 



and 



OT Ap 2 8 2 
dp 



3tt 2 



■P{1 " 



8p 3 

(P-1) 2 
8p 3 



3(1 -x 2 ) In 



1 — Xi 



+ 4(p + l)x? - 6xi 



3(1 - x\) In l±Il + 4( p _ i) x 3 + 6a .l } = ( 34 ) 



1 - x 2 



If we subsitute into (30) the solution to the gap equation (33), we find 



r = %{5* + - 1a 2 

TT Z 3 2 



-l(p+lfxl + (p-lfxl]6 2 }. 



Furthermore, the second order derivatives provide us with the expressions 



d 2 Y 



+ 



(<9A 2 ) 2 27rV 2 L(p+l)(l+x 1 ) (p-l)(l+x 2 ) 



(35) 



(36) 



dp 2 



and 



Ap 2 5 2 



3tt 2 



1 + 



3(P+1) S 
4p 3 



1+Xi 



l-x 2 )ln — 



+ 



-Xi 



P° 
d 2 T 



dA 2 dp 2ir 2 p 2 



In 



2p 3 
1 — Xi 



Xi 
Xi + 



3(p-_l)! (1 _ 4)ln i±- 



4p 3 



x 2 



(p + n\., 3(p + i) ( y + p + i) i , (p-i) 3 ^ 3 ^- 1 )('° 2 -' + 1 )^7) 



2p 3 



In - — — + 2(p + l)xi + 2(p - l)x 2 



(38) 



1 + x 1 1 + x 2 

Our expression of the free energy function eq. (35) coincides with the expression 
derived in [2] if we perform the following modifications: replace the Fermi velocity in 
their formula with one and multiply it by a factor of four. This factor accounts for 
the four pairing configurations, (r,u R ) - (g,d R ), {r,u L ) - (g,d L ), {r,d R ) - {g,u R ) and 
(r, di) — {g, Ul), per momentum for quarks versus one pairing configuration, spin up-spin 
down, per momentum for electrons, where the subscripts R and L refer to the helicities 
of the pairing quarks. 
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3 The Meissner Masses 



In this section we shall consider the coupling of the quark fields to external gauge fields. 
The mean field NJL Lagrangian is modified 

A 2 d 

+ A(-xc75A 2 t 2 X + Xl5^2T 2 Xc)- (39) 

A 1 ^ and represent the gluon and photon gauge potentials respectively. The Lagrangian 
is invariant under SU(3) C x U(l) em transformations. The SU(3) C generators and the 
U(l) em charge operator can be written as 

Ti = \\ h Q = \ + \t,. (40) 
The squared Meissner masses are the eigenvalues of the Meissner tensors, defined by 

(M 2 )g=hmng(X), (41) 

where IT0 (K) is the gluon/photon magnetic polarization tensor given by 

= \ J Q (kBT) 1 dx A J d 3 re- iK * < j[(f, r) j}(0, 0) > c 

= - k -fY, j ^TrViS{P)V)S{P + Kl (42) 

where < ... > c stands for the connected Green's functions, i.e. < AB >c=< AB > — < 
A >< B > and K = (kz,—uj) is the Euclidean four momentum either of the gluon or 
of the photon. The current operator reads J 1 = l^f'^, where the quark-gluon/photon 
vertex with respect to the NG basis (15) takes the form 

r" =( ~ 1 (43) 

where / = 1,2,..., 9 and the tilde denotes transpose. Here we have introduced a new 
set of generators, 77, that mixes the gluons with photon, in accordance with the residual 
symmetry in the presence of the condensate (4) [18], i.e. 

%, — Ti (44) 

for I = 1,...,7, 

T 8 = T 8 cos# + -Qsintf (45) 
9 

and 

T 9 = - -T 8 sin 9 + Q cos (46) 

e 
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with cos 9 = and sin# = — . The condensate is invariant under the residual 

gauge transformations SU(2) x ZY(1), where 577(2) is generated by Ti, T 2 , T 3 and W(l) by 
T 9 . We shall refer to the photon-gluon mixture associated to T 8 as the 8th gluon. 

Each component of the Meissner tensor is a matrix with respect to the generators of 
the gauge group, eqs.(44), (45) and (46). A simplification follows from the observation 
that the isospin part of the quark-photon vertex is the unit matrix with respect to the NG 
blocks ( —r 2 r 3 r 2 = t 3 ) and consequently it contributes to the integrand of the Meissner 
tensor a total derivative, i.e. 



/ 



rf¥ Tr l3 r 3 S(P)MS(P) = 1 1 ^l-TrntMSiP), (47) 



(2tt) 3 u ° y J y ' J (2nfdp 



where M. is an arbitrary matrix. Therefore the isospin term of the charge operator (40) 
can be omitted and we have effectively Q — \ [15, 16]. 

Next we shall prove that the Meissner tensor is diagonal with respect to the indices / 
and I', more specifically that it has the following form 

(MX = (M 2 )J = (M 2 )f = (M 2 )g> = 0, (48) 



and 



(MX = (M^f = (M^f = (MX - (m%, (49) 

(MX = (m'X- (50) 

The relations above follow in a straightforward manner from symmetry arguments, 
in particular the residual SU(2) symmetry. We proceed to divide the set of the nine 
generators into three subsets. Subset I includes the unbroken generators 7i, T 2 and 7jj, 
subset II includes the broken generators T4, %, % and T 7 while subset III includes the 
broken % and the unbroken % generators. Under an SU(2) transformation, % — > u%u\ 
the generators in the subset I transform as the adjoint representation, the generators 
in the subset II as the fundamental one while the generators in the subset III remain 
invariant. Therefore the indices V and I which correspond to generators from different 
subsets do not mix in the expression for (M 2 ) l -J. The mass matrix vanishes within the 
subset I because of the unbroken gauge symmetry. As for generators of the subset II, we 
proceed to relabel its members according to Ji = T 4 , J\ = T 5 , J 2 = T 6 and J 2 = T 7 . 
Introducing 

h=(l) <h = ( ? V (51) 



we have 




(j) a 

€ 



(52) 



and similar expressions for with g = —i(p a , where the block decomposition is with 
respect to color indices. Since J a is symmetric and is antisymmetric, there is no 
mixing between J a and following from the properties of the matrix C = ^7274, i.e. 
C^^C- 1 = -7„ and 

p 3 CS(P)(p 3 C)- 1 = -S(P), (53) 
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a consequence of the symmetry under a combined transformation of a space inversion, a 
time reversal operation and a 75 multiplication (see appendix C for details). The identity 
of the matrix elements (M 2 ) a/3 and (M 2 ) a/3 follows from the equivalence of the doublet 
representation supported by T a and the anti-doublet one by Finally in the subset 
III, the color matrices never mix the red and green components with the blue ones. The 
latter do not carry a condensate. It follows then that we need only to work within the 
red-green subspace, in which % does not contribute and % may be replaced by a factor 

The above argument on the group theoretic structure of the Meissner tensor leads also 
to its relation with the momentum susceptibility introduced in Ref. [15], 



M« = I^ 2 4)aS-' (54) 

Because of the single plane wave ansatz, eq. (4), each of the nontrivial tensors (m 2 )ij 
and (m' 2 )ij can further be decomposed into its transverse and logitudinal components, 
i.e. 

(m 2 )y = A(Sij - q4j) + Bq^ (55) 
(m'% = C(5ij - qiqj) + Dq^. (56) 



In particular, eq. (54) implies that 

~19nA2V + Rn^J F)n ~ f\^{ + Zal) Pin* ^ > 



12p5 2 V 3g 2 > dp Q5 2 ^ 3g 2 ^dp 2 
and 



125 2 V 3g 2 ^dp 2 

It follows from eqs. (57) and (34) at p 7^ that the Meissner tensor of the 8th 
gluon is longitudinal at the LOFF equilibrium. So is, the Meissner tensor of an electronic 
LOFF superconductor. This conclusion of ours disagrees with that of [2], where a nonzero 
transverse Meissner mass was reported. The above property does not contradict with the 
5ij structure of the Meissner tensor of a homogeneous superconductor since it corresponds 
to the trivial solution, p = of eq. (34), where J^j 7^ 0. 

In order to calculate the quantities A, B, C and D we follow the same procedure that 
we used to calculate the free energy in the previous section. In order to improve the UV 
convergence, we employ the technique of [19] and subtract from the polarization tensor of 
the superconducting phase, (42), the correponding expression of the normal phase with 
the same p and 5, which is zero at K = 0, and write 

ng(0) = E / ^Tr^SiP^SiP) - r|'5(P)r;.5(P)], (59) 

where S(P) = S(P) |a=o- Then the weak coupling approximation of the propagator, 
eq.(21), is inserted. We take the z-axis to be parallel to q and carry out the integration 
over the azimuthal angle and the magnitude of p. We are then left with the integration 
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over the meridian angle and the summation over the Matsubara energies. Finally, we 
project out the transverse and longitudinal components and obtain the expressions 

3 n 



A 
B 



d cosO sin 2 0F(0, A), 



dcostfcos 2 9F(9, A), 



C = | ^ d cos sin 2 9G(9, A) , 



and 



where 



and 



D — - 



d cos 9 cos 2 9G(9,A), 



F(9,A) 



m ^ z(v,9,A)iv + 5 e + z(v,9,A) 



G(9,A) = 



3tt 

2g 2 jj 2 A 2 k B T 



( 1 + ^ lm ^7H^Ay 



9vr V 3^ 

At T = 0, the summation over v becomes an integral and we find 

5 2 



F(9,A) 



3tt 2 



51 



2 -f 2+ 9(\S e \-A) 



and 



9tt 2 



i - 



A 2 ^ 



A 2 



A 2 



(60) 

(61) 
(62) 

(63) 

(64) 

(65) 

(66) 
(67) 



We recognize that eq. (66) and (67) are precisely the expressions that Huang and 
Shovkovy obtained for the Meissner masses [13] when their parameter 5/i is replaced by 
an angle- dependent one, \5g\. By carrying out the angular integration as it is described 
in detail in the appendix B, we find 



+ 



6vr 2 
5 2 



25 



at( 1 + ^ 



3A 2 
16p 3 5 2 



In 



1 + 
1 — X\ 



lni±^) 



(68) 



■(P + 1) 5 ^ 



<5 2 



:(p+ 1)^(5x1 -3) 



5p 3 A 2Vr y 1 8p 3 A 2 
- ^(p-l) 5 ^ + ^(p-l) 4 ^(5x 2 -3) 



5p 3 A 2 



8p 3 A 2 



5 



+ 
+ 



1 - 



2^ 
A 2 



67T 2 
J 2 

pA 2 



ft . ° 2\ . 3A " /, 1 i 1 "^ x '2\ 

( 1 + 5 p ) + 8^( ln r^- ln r^) 



8p : 



1 +Xi 

£ 2 

4 p 3 A 2 

<5 2 

4p 3 A 2 



1 + 2V 

(p+l) 4 xi(5x 2 -3) 



(p - lfx 2 {hx 2 2 - 3) 



(69) 
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c = 



and 



9tt 2 



(1 + ^){1 - ^±ll![3(l-,Dln^ + 4(p + l),3_ 6xi 
(P-1) 2 



8p 3 



-3(1 -x z 2 ) In 



l + x 2 
1 - x 2 



+ 4(p - l)x 6 2 + 6x 2 } (70) 



D 



9tt 2 



(1 + ^){1 + ^[3(p + l)(l-^)ln 



P + lr 



+ 



4p 



1 + 
1 — X\ 



+ 4(p+l) 2 ^-6(p 2 + p + l)x 1 



p-1 

4p 3 



-3(p-l)(l-^)ln T 



1 + x 2 



^2 



+ 4(p - l) 2 x'l - Q(p 2 - p + l)x 2 }. 



(71) 



Comparing (70) and (71) with (34) and (37), we confirm explicitly the relation between 
the Meissner tensor of the 8th gluon and the momentum susceptibility, eq. (54). 



4 The Meissner Effect for a Small Gap Parameter 

In the last two sections we derived the expressions for the free energy and the Meissner 
masses of the gluons for a 2SC-LOFF of arbitrary A, 5 and p. In this section we shall 
proceed to address the issue of the chromomagnetic instability at equilibrium. Initially we 
shall consider the case A << 5, -the gap parameter much smaller than the displacement. 
This is the situation near a second order phase transition to the normal phase. For a 
LOFF state consisting of a single plane wave at T = this occurs at [2] 

5 = 5 C ~ 0.754A , (72) 

which is also the maximum for the LOFF condensate, eq. (4). Although it has been 
argued that the single plane wave ansatz is not energetically favorable compared with 
the multi plane wave ansatz in this neighbourhood and that the transition to the normal 
phase may be a first order one [21], we shall still pursue the small gap expansion since 
it provides a full analytic treatment and consequently can be very instructive. We shall 
focus on the LOFF state at T = 0. 

Let's expand the free energy, eq. (30) to quartic power of the small parameter, A = 
A/5, 



r = 2p 2 5 2 

with, 



7T 2 



a A 2 + 1(3 A 4 } (73) 



and 



a = ln^-l + ^ln(p+l) + ^-iln|p-l| (74) 
A 2p 2p 

fi = W=ry (75) 
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The corresponding expansions of the quantities A and B read 



A 



B 



6tt 2 

g 2 p 

' 6tt 2 



(a A 2 + aiA 4 + .. 



(6 A 2 + 61A 4 + 



where 



a 



Oi = 



8p ; 

16p 3 ^p — 1 p + 



P + 1 
p-1 



t)- 



1 



32p 3 L(p-l) 2 (p+l)< 



6o = 4^( 



3 ,2-p 



4p2 Vp2 



1 P P 



and 

Similarly, we find 



+ ■ 



P + 



t) 



+ 



+ ■ 



16p 3 L(p-l) 2 (p+1) 2 J 16pL(p+l) 3 (P-1) 3 



2^2 



c = 



D = 



9 P 
9tt 2 

2—2 

9tt 2 



(l + ^)(coA 2 + Cl A 4 + ...), 

,2 



(l + ±j) (doA 2 + d x A 4 + ...). 



where 



c = 



Cl 



32p 3 
d 



4p 3 
2 



(2p - In 



P + 1 



2a , 



p+1 p-1 

3 /2-p 
2p 2 V _ 



(p 
P- 



il„: 

P P 



I) 2 (P-1) 2 



and 



+ 



3 



16p 3 Lp+l p-1 (p + 1) 2 (p-1) 2 (p+1) 3 (p-lf 



(76) 
(77) 

(78) 
(79) 

(80) 
(81) 

(82) 
(83) 

(84) 
(85) 
(86) 

(87) 



The validity of the relations Cq = 2ao and d = 26 is not restricted to zero temperature 
and the weak coupling approximation employed here. If we expand the integrand of the 
Meissner mass tensor (59) to the order A 2 , we find that the loop integrals for ao and bo 
are related to those for c and d through integration by parts [15]. The same argument 
also leads to identical criteria for the chromomagnetic instability of the 4- 7th gluons and 
the 8th gluon near T c for g2SC. Minimization of the free energy with respect to A and p 
yields 

A 2 = 4(p 2 -l)(l-A) ; ( 88 ) 
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and 

P = p c + \-^-A 2 , (89) 
where p c ~ 1.20 is the solution of the equation [2] 

J- In ^±1 = 1 (90) 
2p c Pc-1 

and 

5 1 

= — . ~ 0.754. (91) 

Ao 2^1 

The free energy at the minimum reads 

4p 2 



ruin o 
7T 



(p 2 -l)(5-<y 2 (92) 



and both branches of eq. (23) are gapless here. The solution satisfies the stability condi- 
tion (13). Substituting it into eqs.(76), (77), (82) and (83), we obtain 

a 2 u 2 A 4 a 2 u 2 A 2 

A = 9 9 - > 0, B — 9 ^ 9 > (93) 

96tt 2 (p 2 - l) 2 - ' 8tt 2 p 2 c - 1 - K ' 

and 



C = 0, £ > = ^-(l + ^-)_>0 (94) 



6tt 2 V 3^ 2 ^p 2 -l 

Therefore the chromomagnetic instabilities are removed completely for small gap param- 
eters. 

The situation for a multi-plane wave ansatz for the LOFF state will be better if the 
small gap expansion remains numerically approximate. As we have seen for the single 
plane wave ansatz, the Meissner tensors are purely longitudinal to A 2 and the logitudinal 
components are positive at the equilibrium. It follows that if the LOFF state (3) contains 
three linearly independent momenta, the Meissner tensor will be positive definite, since 
it is additive with respect to different terms of eq. (3) to order A 2 . 



5 Numerical Results 

In the last section, we have shown analytically that the chromomagnetic instability is 
absent in the region near the outer edge of the LOFF window at zero temperature. In 
this section, we shall explore the interior region in order to determine if the absence of 
the chromomagnetic instability can be extended to the entire LOFF region at T = 0. 

Before presenting our numerical results, we shall emphasize some general features of 
the solutions of eqs. (33) and (34). A trivial solution to eq. (34), p = 0, corresponds to a 
homogeneous 2SC and there is always a solution to eq.(33) with A = A for < 5 < A . 
The expression for the free energy then becomes 
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Since the excitation spectrum, eq. (23), is fully gapped in this case, we shall refer to the 
correponding state as the BCS state even in the presence of a displacement, following the 
convention of Ref. [2]. We observe that the BCS state becomes meta-stable relative to 
the normal phase for 

5 > ^=A . (96) 

For the displacement within the interval ^ < 5 < A , there is another solution at p = 
with 

A = y/A (26- A ) (97) 
which leads to the following expression for the free energy 

r = 4(25-A ) 2 >0. (98) 

This solution corresponds to the Sarma state that is unstable in the grand canonical 
ensemble for a given 8. In a canonical ensemble, it may be the only charge neutral state 
for certain pairing strength [7] but suffers from chromomagnetic instabilities [13]. 

In what follows we shall restrict the terminology 'LOFF state' to the nontrivial solution 
to eq. (34) that satisfies 

(P + 1)2 [3(1 - x\) In + 4(p + l)x\ - Q Xl 

1 1 — X\ 



3 



8p- 



-3(1 -a; 2 ) In — — - + 4(p - l)x 6 2 + 6x 2 = (99) 

Xi 



8p 3 

together with eq. (33). It follows from eq.(33) that 

A < A (100) 

at the solution and eq.(99) implies that 

y<l + P- (101) 

Otherwise, X\ = x<i = and eq.(99) cannot hold. Physically, eq.(100) is the consequence of 
the reduced phase space available for pairing and eq.(lOl) is the condition for the existence 
of gapless excitations, which is necessary to make the net baryon number current vanish. 
Combining eq. (100) and (101), we find 5 < A for the LOFF state. 

For a small gap parameter, an analytical solution was obtained in the last section, 
which corresponds to a stable state with lower energy compared with the BCS and the 
normal states and is free of chromomagnetic instabilities. For a general gap parameter, 
however, eqs. (33) and (99) can only be solved numerically. Our strategy is to find 
solutions of (99) for A/ 5 for a given p and subsequently to calculate the corresponding 
S/Aq in terms of (33). Concurrently, we examine the local stability according to eqs. (13) 
and the competition between the LOFF state we find and the BCS state of the same 
5/Aq. The Meissner tensors are calculated as well. We observe that the quantities p, 
5/A , the free energy difference, the Meissner masses and the stability determinant of eq. 
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13), at the LOFF solution are single- valued functions of the normalized gap parameter 
A/A . These functions are plotted in Figs. 1-7. 

The diquark momentum of a LOFF state, normalized by the displacement parameter 
versus the normalized gap parameter is shown in Fig.l and the corresponding ratio of the 
displacement parameter to the BCS gap is shown in Fig.2. In Fig.3, we display the free 
energy difference between a LOFF state and a BCS state at the same 8, normalized by 
the magnitude of the BCS energy (95) at 5 = 0. Combining the three figures, we find the 
LOFF window, in which the LOFF state is of the lowest energy, 

max; 

(102) 

where <5 max = 5 C ~ 0.754A and 5 min ~ 0.706A , slightly below the threshold (96). The 
corresponding momentum of the diquark increases monotonically from p ~ 1.20 at 5 max 
to p ~ 1.28 at 5 m i n . Our results confirm those of Ref. [2]. In addition, we found that the 
normalized gap parameter, A/A varies monotonically from at 5 max to 0.242 at <5 mm . 
So the LOFF window on all figures corresponds to the domain < A/A < 0.242 of the 
abscissa. 

The numerical results of the Meissner tensors, Figs. 4-6 are most interesting. The 
numbers shown are normalized by the Meissner mass square of the BCS state at 5 = 0, 
i.e. 

(m% = (103) 

and 

(m% = D Q 8ij (104) 

with 

All squared Meissner tensors are nonnegative within the LOFF window and the region 
free from the chromomagnetic instabilities extends beyond it. While the Meissner tensor 
associated to the eighth gluon is aways logitudinal, Fig. 6, the transverse component of 
the Meissner tensor associated to the 4- 7th gluons, Fig. 4 is aways much smaller than its 
logitudinal component, Fig. 5. The stability of the solution against virtual displacements 
of A and p within the LOFF window is guaranteeted by the positivity of D, plotted in 
Fig. 6 and the positivity of the determinant in eq.(13), plotted in Fig. 7. The non- 
smooth behavior of the function D and the determinant of Figs. 6 and 7 correspond to 
the transition of the LOFF state from the region where both branches of (23) are gapless, 
x 2 > 0, to the region where only the Ez branch is gapless, x 2 = 0. The derivative of 
x 2 with respect to both A and p become singular there. When the gap parameter is 
expressed in terms of p and x 2 > according to eq. (32) and the function D and the 
determinant are expanded according to the powers of x 2 , the linear term is sensitive to 
the singularity and gives rise to the non-smooth behavior. On the other hand, the same 
type of expansion of the functions A and B does not contain a linear term in x 2 and thus 
the transition behavior of A and B is much more smooth. 

As the LOFF momentum goes to zero, p — > 0, we have A — > Ao and 5 — > Ao, a result 
can be established analytically. Therefore the LOFF state, the BCS state and the Sarma 
state joins at the point ( A , A ) in the three dimensional parameter space of A, S and 
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p. We also observe from Figs. 4-7 that the chromomagnetic instability and the instability 
against a small variation of parameters show up for sufficiently low values of p. Therefore 
the instability free region of a LOFF state is disconnected from the g2SC state in the 
parameter space. 

6 Concluding Remarks 

In this paper, we have systematically examined the Meissner effect in a two-flavor color 
superconductor with a LOFF pairing. We found that within the LOFF window at T = 0, 
that is within the range of the displacement parameter where the LOFF state is ener- 
getically favored, all the squares of the Meissner masses are non-negative. Although the 
positivity of the longitudinal component of the Meissner mass tensor associated with the 
eighth gluon follows from the stability of the LOFF state against a variation of the LOFF 
momentum, p, the positivity of the Meissner tensor for the 4- 7th gluons appears acci- 
dental regarding the pairing ansatz (4). We have also found that the Meissner tensor of 
the 8th gluon is always longitudinal. This is also the case with the Meissner effect in an 
electronic LOFF superconductor. 

The charge neutrality condition has to be implemented before we can claim that the 
LOFF state is the ground state of quark matter at moderately high baryon density. In 
the previous letter [15], we showed that the chromomagnetic instability of the 8th gluon 
pertaining to g2SC implies a state of lower free energy with a nonzero diquark momentum 
and without offsetting the charge neutrality. This is not sufficient since the LOFF window, 
with the diquark momentum comparable to the Fermi momentum displacement, is quite 
far away from the g2SC and is very narrow in the three dimensional parameter space. It 
remains to be seen if quark matter within the LOFF window can be made neutral. In 
a hypothetical world where the sum of the electric charges of the u and d quarks is a 
small number, we can find within the framework of the weak coupling approximation an 
analytic expression for an electrically neutral LOFF phase that is located near the upper 
edge of the LOFF window and is energetically favored over both the neutral normal phase 
and the neutral BCS phase. Also it is free of chromomagnetic instabilities. The details 
will be presented in a forthcoming publication [20]. Notice that the neutral phases being 
compared do not have the same displacement parameter 5 since the charge neutrality 
constraint adjusts S differently for different phases. In the real world, the sum of the 
charges of the u and d quarks is 1/3. The weak coupling approximation may be marginal 
and one may have to seek numerical solutions of eq. (33), eq.(96) and the charge neutrality 
constraint simultaneously. 

We also notice that, if the pairing force is mediated by the long range one-gluon 
exchange of QCD instead of the point coupling of NJL model, the LOFF window will 
be at least five times wider under the same Ao {S max ~ 0.968Ao and 5 m i n < 0.706Ao)[6], 
making the implementation of the charge neutrality more likely. 

In this paper we have considered only the simplest LOFF ansatz that consists of one 
plane wave and pairs u and d quarks. It was suggested that multi plane wave condensates 
are favored more near the upper edge of the LOFF window [21]. Numerical analysis also 
suggests that the s quark cannot be ignored in pairing at T = 0. A more complicated 
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structure of the condensate is required for a realistic LOFF pairing in quark matter and 
may not be accessible via analytical means. In any case, the absence of chromomagnetic 
instabilities for the simple 2SC-L0FF is encouraging and promotes the LOFF state as a 
potential candidate of the ground state of quark matter at moderately baryon densities. 
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Appendix A 

In this appendix, we shall first derive the exact expression of the quark propagator in 
the superconducting phase for 5' = 0. Subsequently we shall perform the weak coupling 
expansion. Finally, we shall approximate the determinant of the quark propagator which 
enters the free energy function (20) in a similar manner. 

As r 3 and A2 are the only isospin and color matrices appearing in the inverse quark 
propagator (18) at 8' = 0, they can be treated as a c-numbers. Taking the square of the 
inverse propagator, we have 

[S-\P)f = (tu-5r 3 ) 2 + p 2 -p 2 -q 2 -A 2 X 2 2 (106) 
+ 2p 3 [//(w - 5t 3 ) - p-q\+ 2Api7 5 (-ry • q + £74) A 2 . 

In terms of the matrix 

N(P) = ( tl y-5r 3 ) 2 + l 2 2 -p 2 -q 2 -A 2 Xl 

- 2p 3 [ / u(iz/-5r 3 ) -p-g] -2Ap 1 7 5 (-i7-g + / u7 4 )A 2 , (107) 

we obtain the exact expression of the quark propagator 

«p)=qgp, (108 ) 

where 

D(P) = [{iv-5r 3 ) 2 + p 2 -p 2 -q 2 -A 2 \ 2 2 ] 2 -A[p{iv-5T 3 )-p-qf + 4A 2 (fi 2 -q 2 )\ 2 2 , (109) 

and \ 2 projects into the subspace of the red and green colors. 

In order to perform the weak coupling approximation, we regard p — p,, v, 5 and q as 
small quantities of comparable magnitudes and maintain only the leading order terms in 
the numerator and the denominator of (108). We find that 

D(P) ~ Afi 2 [-(iu -Sr 3 -p- q) 2 + (p- pf + A 2 \ 2 2 ] (110) 

and 

S-\P)N{P) ~ -2p 2 H7 • p + p 374 )[|p - p 3 <?1 -Ji + psiiu- 5r 3 ) + AA 27 57 4 ] (HI) 
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Upon substitution of (110) and (111) into (108), we derive the approximate quark prop- 
agator expression (21) that appeared in section II. 

The determinant of S~ 1 (P) may be calculated similarly. Let's write 

detS-^P) = det^S-^P)] 2 ] ~ (2p) 24 det^(P), (112) 

where 

d(P) = p-p + p-i{iv -8r 3 -p-q) + AA 2 pi7574- (113) 

Now by decomposing the 48x48 matrix d(P), according to the eigenvalues of r 3 and A 2 , 
i.e. 

detd(P) = 1] detd(P) | T3 =±i,a 2 =±i,o (H4) 

we have 

detrf(P) | r3 =±i,A 2 =i= [(p -fL) 2 + A 2 -(iv T 5-p- q) 2 ]\ (115) 
detd(P) U =±1) a 2= _i= [(p - p) 2 + A 2 - {iv T 5 - p ■ q) 2 ] 4 (116) 

and 

detd(P) |. 3 =±i,A 2 =o= [{p - pf - {%v T S - p ■ q) 2 f. (117) 
Eq. (22) is then established. 

Appendix B 

The integrand of the angular integrals for the free energy function (28) and that for the 
Meissner tensor, eqs. (66) and (67), contain terms with 9{\5q\ — A). It confines the integral 
within the region in the phase space where either or both branches of the excitation 
spectrum, Ep of (23), become gapless. This occurs when A < (1 + p)8. Otherwise the 
integration domain specified by 9(\5$\ — A) is null. There are three different cases to be 
considered: 

1) |1 - p|5 < A < (1 + p)8: 

Only the branch may be gapless and we have 

y dcos00(|£ fl |-A)(...) = J ' " dcos0(...). (118) 

2) p < 1 and A < (1 - p)5: 

Again we may have only one gapless branch, E~, but instead 

J 1 d cos 66(\6 e \ - A)(...) = J* dcos9(...). (119) 



because > I. 

3) p > 1 and (p - 1)5 < A < (p + 1)5: 

Both branches, E~ may be gapless and we have 

. S-A 



J dcos99(\6 e \ - A)(...) = J pS dcos9(...)+ fdcos9(...). (120) 

The indefinite integrals of the integrands of the free energy function and the Meissner 
tensors are all elementary functions. 
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Appendix C 

In this appendix, we shall supply the details supporting the statement in section III that 
the identity (53) for the quark propagator follows from the invariance of the 2SC-LOFF 
under a combined transformation of a space inversion, a time reversal and a 75 multipli- 
cation. We shall work with the canonical formulation at T = though the generalization 
to T 7^ is straightforward. 

The Hamiltonian operator corresponding to the mean field NJL action (9) reads 



A 2 r 

H MF = Vl— — + / d 3 f[^7 • - ^74/^ 

A(-^75A 2 r 2 ^e- 2 ^ + ^75A 2 r 2 ^ c e 2 ^)]- (121) 

Introduce an anti-unitary operator U which implements a space inversion, a time reversal 
and a 75 multiplication. The transformation laws of the field operators in Heisenberg 
representation are given by 

U^{r,t)U~ l = i^i-r, -t) 

Ui)\r,t)W x = -i^(-r,-t)i2. (122) 
The corresponding transformations of the NG spinors are 

W*(f,*)W -1 = wy2*(-r,-*) 

U^if^U- 1 = -i^(-r,-t) l2 p 3 , (123) 

following (5) and (15). The phase factor associated to the time reversal has been fixed in 
accordance with phase convention of the condensate such that 

UHmfU- 1 = H MF - (124) 

For a state | > that is invariant under U, say the ground state and an arbitrary operator 
O, we have [22] 

< \0 ] \ >=< \UOW x \ > . (125) 

The Hamiltonian (121) is also invariant under a unitary transformation, U, that generates 
a translation and a rephasing according to 

U^(r,t)U* = e* s il>*(r + a,t). (126) 
Consider the quark propagator in coordinate space, 

S aP (f,t) =< |^(f,t)^(0,0)| > 6{t)- < 1^(0,0)^(^)1 > 0(~t). (127) 
It follows from (123), (125), (126) and time translation invariance that 

Safi(r,t) = -(p 3 C) a ,S a ,(r,t)(Cp 3 )p a , (128) 
which, upon a Fourier transformation, gives rise to the identity (53). 
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Figure 1: The dimensionless diquark momentum p of the LOFF state as a function of the 
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Figure 2: The normalized displacement parameter, 5/Aq, of a LOFF state as a function 
of the normalized gap parameter A/A . 
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Figure 3: The free energy difference between a LOFF state and a BCS state as a function 
of the normalized gap parameter A/A . The coordinate is normalized to the magnitude 
of the free energy of the BCS state, eq.(94), at 5 = and the LOFF state is favored when 
the value is negative. 
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Figure 4: The normalized transverse Meissner mass square of the 4-7 gluons, A/A , as a 
function of the normalized gap parameter A/A of a LOFF state. 
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Figure 5: The normalized logitudinal Meissner mass square of the 4- 7th gluons, B/A , as 
a function of the normalized gap parameter A/A of a LOFF state. 
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Figure 6: The normalized logitudinal Meissner mass square of the 8th gluon, D/D , as a 
function of the normalized gap parameter A/A of a LOFF state. 




Figure 7: The stability determinant in eq. (13) (arbitrary unit) of a LOFF state as a 
function of the normalized gap parameter A/A . 
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